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I. INTRODUCTION 



Strong quantum electrodynamics (QED) significantly differs from weak perturbative 
QED. Pairs of bosons and fermions are produced from external electric fields exceeding 
the critical strength Ec — m^c^/fiq. Vacuum fluctuations polarize the vacuum and the in- 
stability therefrom leads to pairs of particles [1, 2]. This pair production of charged particles 
and antiparticles is a nonpertubativc effect of QED. The recent development of X-ray free 
electron laser technology may make a test possible for pair production and attracts interests 
in this strong QED [3, 4]. Other natural phenomena associated with strong electromg- 
netic fields are astrophysical objects. For instance, pulsars have magnetic fields comparable 
with or above the critical strength, and the magnetic fields of magnetars exceed the critical 
strength by several order of magnitude. Likewise, charged black holes, especially the ex- 
tremal ones, can have electric fields far greater than the critical strength. Recently, Ruffini 
et al. [5, 6, 7] proposed a mechanism for gamma rays burst through pair production and 
annihilation from charged black holes. 

In this paper we introduce the recent canonical method for pair production [8, 9], which 
can be applied to inhomogeneous electromagnetic fields such as charged black holes. In 
fact, the electric fields by charged black holes are inhomogeneous in that the direction is 
radial and the field strength depends on the radial distance. The pair-production rate by the 
proper time method was derived for static uniform fields [1, 2]. However, this formula was 
used to calculate the pair- production rate by charged back holes [5, 6, 7, 10]. Thus a field 
theoretical method is needed to apply even to such inhomogeneous fields. Our canonical 
method is an extension of the early idea that particle production is related with tunneling 
processes in the space-dependent gauge [11, 12, 13, 14, 15]. In the space-dependent gauge 
of a static uniform electric field, each mode of the Klein-Gordon or Dirac or vector field 
equations has a potential barrier. It is the observation that the tunneling probability has 
relation with pair production and the no-tunneling probability is related with the vacuum 
persistence. Further, the tunneling probability is determined by the instanton action includ- 
ing all quantum corrections. Thus this canonical method can readily be applied to static 
inhomogeneous electric fields, in particular, charged black holes. 

The organization of this paper is as follows. In Sec. II we review the vacuum polarization 
and pair production in proper time method. In Sec. Ill we use the time- and space- 
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dependent gauges to obtain the pair-production rate by a uniform electric field. We compare 
the similarity and differences between two gauges. It is shown that particle production is 
related with the tunneling probability and the vacuum persistence is determined by the no- 
tunneling probability. In Sec. IV we discuss the spin effect on pair production and clarify 
the terminology of pair-production probability, pair-production rate and the mean number 
of pairs. In particular, the tunneling probability is given by the instanton action including 
all quantum corrections. In Sec. V we study the Klein- Gordon equation first in Minkowski 
spacetime and then in charged black hole background. It is shown that there are quasi- 
stationary states describing purely outgoing waves with complex frequencies corresponding 
to pair-production states. Finally, wc apply the canonical method to calculate the pair- 
production rates by charged black holes. 



II. VACUUM POLARIZATION AND PAIR PRODUCTION 



At the dawning of quantum field theory, Heisenberg and Euler calculated perturbatively 
the effective action for uniform electromagnetic fields [1]. Later, Schwinger found the one- 
loop effective action for fermions of spin 1/2 in uniform electromagnetic fields by using the 
proper time method [2]: 
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The one-loop effective action can be written in the form 
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which has simple poles at s„ = (n7r)/(eXj), (n = 1,2,---). These poles contribute the 
imaginary part 

2I-£.-J^i:^e-»'-.cth(§„.). (5) 
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The real part of the one-loop effective action gives rise to the vacuum polarization 

x(esXi)2'=-^coth(esX^)} - 1 - -{esfj^ 

--^+J^ (t^)'^^^ - ^"f + 7(E • B)^] + • • • (6) 
45 HTrhc' mc^ ^ ^ M ^ ^ 

For a pure electric field (B — O^Xr — 0,Xj = E), Eq. (5) leads to the pair-production 

rate 

— e"^^. (7) 
In the case of bosons, the pair-production rate becomes 

2Im£^ = ^ E(-l)"^^(-) e-^. (8) 

These imaginary parts determine the vacuum persistence (vacuum-to- vacuum transition) 

|(0,out|0,in)|2 = e-2^^i°^^% (9) 

and thus lead to pair production of particle and antiparticle. It should be remarked that 
the pair-production rates (5), (7), (8) and the vacuum polarization (6) are valid only for 
uniform fields. These formulae may be used as long as fields vary slowly in region of interest. 
There are many situations, such as charged black holes, where the fields are strong but not 
uniform. In the next sections, using a canonical method, we shall derive the pair-production 
rates valid even for inhomogeneous fields. 



III. CANONICAL METHOD FOR PAIR PRODUCTION 

Pair production of charged particles and antiparticles can also be understood in canonical 
theory. An advantage of the canonical method is that it has a direct meaning of pair 
production compared with the path integral methods. In the canonical method one directly 
solves either the Klein-Gordon or Dirac equation. An interesting point is that one can easily 
extend the result for bosons to fermions by taking into account the Pauh blocking effect due 
to spins. The Klein-Gordon equation for bosons is given by 
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In this section we focus on the uniform electric field, which is determined by either vector 
potential or scalar potential 

E=-|a-VAo. (11) 
A uniform electric field has a time-dependent vector potential 

A = -Etk (12) 

or space-dependent scalar potential 

^0 = -Ez. (13) 



A. Time-dependent Gauge 



In the time-dependent gauge (12) with the electric field in the 2;-direction, the Klein- 
Gordon equation for spin-0 bosons now takes the form 



9t -dl- {d, + iqEtf + 



$ = 0. 



(14) 



The field is decomposed into Fourier mode, $ = e^'^^ '^^k; and results in a one-dimensional 
equation in the time-direction 



-dl - {K + qEtf - (m^ + ki) 0k = 0. 



(15) 



The above equation describes a particle moving over an upside down oscillator potential 
with positive energy. This equation is an analog of fields in a time-dependent spacetime 
background. We may apply the canonical method by Parker and DeWitt to calculate the 
amount of particle creation [16, 17, 18]. The idea is that an incoming positive frequency 
solution at the past infinity scatters by the potential and separates into an outgoing positive 
and negative frequency solutions at the future infinity. The negative frequency solution 
leads to particle creation. Indeed, Eq. (15) has an incoming positive frequency solution at 
t — > — oo 



0k,m(^) 
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where 
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T = j2qEt, K= Jl + e-^'^^'^x - e^'^^^'x . 



(16) 
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The incoming wave can be written in another form 

1 



0k,m(^) 



\U - -)E{-ax,^,T) + i)£;*(-ak^,T) 



(18) 



In fact, the incoming wave is a superposition of the outgoing positive and negative frequency 
solution 

0k,m(^) = /^k0k,o«t(^) + i^0k,o«t> (19) 

where 0k,out is the outgoing positive frequency solution given by 

1 
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Equation (19) is equivalent to the Bogoliubov transformation between the operators in two 
asymptotic regions 

bk,out = fJ'\J>k,in + ^k^k,o«t) (22) 

The number of created particles in pairs in the future infinity is given by 



(0,m| E^Ut^k,o„*|0,in) = Y: kkr = Ee-'""''^ 

k kx kx 



(23) 



Here the mean number of created pairs in the transverse mode kj^ is Akx = |i^kP < 1- Then 
the vacuum persistence or vacuum-to- vacuum transition is given by 

1 
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k, i/^kl 



(24) 



where 



;,k= (l + e-2™'^x)V2. (25) 
Thus the vacuum persistence is related with the imaginary part of the effective action as 



I (0, out|0, in) |2 = exp [- ^ ln(l + e-^^i'x ) 

kx 



exp 



-2VTlmjCi 



(26) 



Now the imaginary part of the effective action for bosons per volume per time is given by 
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We may interpret Eq. (19) quantum mechanically. In the space with x — —t, (f)]s.,out corre- 
sponds to an incident wave and 0k to a reflected one. Thus the transmission coefficient 
gives the vacuum persistence 

2 1 
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kx) 



1 + e-^'"''^± 

and the mean number of created pairs is given by the reflection coefficient 



t'k = e 



-27rak I 



(28) 



(29) 



B. Space-dependent Gauge 



In the space-dependent gauge, i.e. the Coulomb gauge (13), the Klein-Gordon equation, 
after the mode-decomposition $ = e'^^^-^'^-^~^^^(f)u,kj_, separates into the mode equation 



-dl - (cu + qEzy + (m^ + k^) 



ya;kx 



0. 



(30) 



This equation describes a particle moving under the potential barrier in contrast with the 
over-barrier in the time-dependent gauge. The tunneling wave function is given in terms of 
the complex cylindrical function 



-ujkAO = C^(«kx,C), 



where 

This wave function has the proper asymptotic forms at z —>■ —oo 



and at z — > -|-oo 



Here the coefficients are 



,kx(e^ oo) = Ckxy||e+2«'. 



^kx = V 1 + e'™''^ , 5k^ = ^ce™''^ , Ckx = c. 



Then the tunneling probability is 



1 + e 



27rak| ' 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



and the probability for no-tunneling 



pn6 



Note that the instanton action for tunneling 



1 + 



(37) 



V = ^ / dz^J{w? + ^_)-{u + qEzf = TTOk^ , (38) 

after taking into account multi-instantons and multi-anti-instantons contributions, leads to 
the total tunneling probability 

OO -1 

PL = E(-l)"^^e-"^''x = .^^^^ (39) 

n=l -I- -t- c 

and the total probability for no-tunneling 

OO -I 

n=0 J- + e -L 

The idea of Refs. [8, 9, 11, 12, 13, 14, 15] is that the tunneling probability is related with 
pair production and the no-tunneling probability with the vacuum-to- vacuum transition. As 
in the case of the time-dependent gauge, the vacuum persistence is given by 

|(0,out|0,in)|2 = JJp^l = exp[-2yrim£^]. (41) 

kx 

Finally, the imaginary part of the effective action for bosons per volume per time is given 
by 

21-4 = ^^Eln(l + e->.) = pL,|(-l)"«(^)'--"V=S^. (42) 

In the fermion case, fermion pairs in the same state are prohibited by the Pauli exclusion 
principle, so the tunneling probability is limited to 

= e-^''''^± , (43) 

and the no-tunneling probability to 

P^f = 1 - e-^'^"''x . (44) 
Thus the imaginary part of the effective action for fermions per volume per time is given by 
2Imi:{ = -^Eln(l - e-^'^^) = E(-) ■ (45) 

kx ^ ' n=\ 
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IV. SOME REMARKS ON CANONICAL METHOD 



In this section we discuss some issues of pair production in the canonical method. In 
Sec. Ill we derived the pair-production rate of spin-0 bosons and spin-1/2 fermions. Here 
we introduce a recent result extending the canonical method to spin particles [9]. We 
also clarify the different terminology for particle production in literature. Finally, we show 
that the canonical method based on the instanton action can readily be generalized to 
inhomogeneous fields. 



A. Spin and Scattering Processes 



To treat the bosons with spin a > 1 and fermions with a > 3/2, one has to solve 
the Klein-Gordon or Dirac equation coupled with the vector field. Then the Klein-Gordon 

equation for spin cx-particles takes the form 



After mode-decomposition, it reduces to 



(46) 



-di - (u; + qEzf + {mf + kl) + liaqE 4>^^^,^ = 



(47) 



The tunneling wave function is now given by 



W(0 = c^(«ki,0> 



(48) 



where 



Denoting the asymptotic form of the wave function 
we find the wave function at = — oo 



(49) 



(50) 



(51) 



and at ^ = -|-oo 



0a;kx = C'kxV'il^kx- 



(52) 
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In terms of group velocity, the flux conservation reads for bosons 



and for fermions 



l^kxT + l^kxT — l-^kxT- 



Therefore, the reflection coefficient for bosons is 



-Bkx 


2 

= 1 - 


Ckx 


^kx 




^k, 



and for fermions 



^kx 


2 

= 1 - 


Ckx 


-Bkx 




-Bkx 



(53) 



(54) 



(55) 



(56) 



This result makes the difference of pair production between bosons and fermions. 



B. Pair Production Rate and Mean Number 



Different terminology for pair production has been used in literature. In some cases one 
has to take care of the different quantities. Let us denote 



w — 2Im>Cp 



(57) 



per volume per time, then in the case of both electric and magnetic fields this quantity 
becomes 



Wf 



= ^4^E(-ir^-cosech(!^)e-^, 
2(27r)2 ti ri ^ E ) 



2(27r) 

/ \o > -coth — — )e iB 
(27r)2 ^^n y E ) 



(58) 
(59) 



Schwinger [2] used the term pair-production probability for or wj, and Itzykson and 
Zuber [19] the pair- production rate, whereas Nikishov [11] used only the first term of the 
series to denote the pair-production rate. 

The mean number of created pairs can be derived in the following way. Let Pq denote the 
probability for the vacuum-to- vacuum transition and P\ denote the probability for one-pair 
production, then for the boson case one has the probability conservation 
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+ •••), 



(60) 



and the mean number of created pairs 

= Po {Pi + 2P' + 3Pf + •••). (61) 

Therefore, one finds the relations 

In the case of fermions, one has the probabihty conservation 

l = Po(l + Pi) (64) 
due to the Pauh blocking, and the mean number of created pairs 

Af^ = PoPi. (65) 

Thus one obtains the relations 

Po = 1 - A/} = P"-^, (66) 



C. Instanton Method 

The space-dependent gauge has an advantage of readily being applicable to static inho- 
mogeneous fields. An inhomogeneous electric field, for instance along the z-direction, has 
the time-dependent gauge 

Az = J E{z)dt = E{z)t. (68) 

So it is difficult to apply the canonical method in the time-dependent gauge. Now, in the 
space-dependent Coulomb gauge, the mode-decomposed field equation takes the form 

-dl + Q{z)]<l>iz)^0. (69) 

Hence the problem reduces to a one-dimensional tunneling problem to which we may apply 
the idea of the canonical method. In general, there are two asymptotic regions with Q{z) < 
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and a potential barrier Q{z) > in between asymptotic regions. Then the no-tunnehng 
probabihty for the vacuum persistence for bosons becomes 

and for fermions 

p/ = 1 - (71) 
where S is the total sum of all order instanton contributions 

oo 

S^T.Si2k). (72) 

k=0 

The 0-loop or classical instanton action is 

2S^o) = jQ"\z), (73) 
and 2S{2k) is the /c-loop correction to the instanton action. 

V. PAIR PRODUCTION FROM CHARGED BLACK HOLES 

We now turn to the pair production by charged black holes. As a charged black hole we 
consider the Reissner-Nordstrom black hole with the metric 

ds" = -g{r)de + ^ + r'^d^l (74) 

g{r) 



where 



Oir) - 1 - ^ + (75) 



Here M and Q are the mass and charge of the black hole. The charge Q also generates an 
electric field in the radial direction given by the Coulomb potential 

Ao = ^, Ai^Q. (76) 
r 

Then the Klein-Gordon equation in the black hole background takes the form 
1 r d .qQ)2 1 d ( ^ , .0 , L2 



It is convenient to introduce the tortoise coordinate 

dr , , , 



[^- = r + - — "^-^ — -ln(r-r+)-- — -ln(r-r_), (78) 

J g[r) (r+— r_) (r+ — r_) 

12 



where r± — M ± \/M'^ — are the outer and inner horizons. The field can be expanded 
into the spherical harmonics 



^(i,x) = e-*^yi„(^,(^), (79) 
r 

where it is assumed Re(a;) > in all cases and Im(a;) = for the bounded states of stable 
motions and Im(a;) > for the unbounded states of unstable motions. Then the radial 
equations for the mode-decomposed field (79) take the form 



+ Vi{t) %{r)^e'^i{T), (80) 



where 

12 



V,(r) = y(r){m' + + ^ - ^) - - 5^)' + (.^ - m% (81) 

and 

Note that Eq. (80) has the same form as the Regge- Wheeler equation for the axial pertur- 
bations of black holes [20, 21]. 



A. Minkowski Spacetime 

In the Minkowski spacetime with M = and Q = in g{r), the effective potential for 
the radial motion reduces to 

l/,W = ^ + ^a±ilf^. (83) 

Note that the effective potential (83) of the Klein-Gorodn equation is the same as that of 
the non-relativistic system with a charge 2ujQ and effective angular momentum /'(/' -|- 1) = 

+ + -q'Q'. (84) 

The characteristics of states depend on the relative magnitude of charge qQ and angular 
momentum Z-l-1/2. The states for \qQ\ < 1 + 1/2 describe stable motions of bound states and 
the hydrogen-hke atoms with Z < (l + l/2)/a belong to this class, where Z is the nucleus 
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charge and a — e^ the fine structure constant [22, 23]. On the other hand, for qQ > 1 + 1/2, 
as in the case of heavy atoms with Z > (l + 1/2) /a, the so-called 'centrifugal' potential 
(83) becomes attractive and the effective angular momentum (84) takes a complex value. 
Moreover, the electric field is sufficiently strong enough to produce pairs of charged particles. 
There are certain quasi-stationary or meta-stable states describing the pair production. 

We consider the interesting case of pair production by the Coulomb field when \qQ\ > 
The centrifugal potential becomes attractive rather than repulsive, and there may exist 
quasi-stationary states with complex frequencies uj 

Re(t^) > 0, Im(a;) > 0. (85) 

Then the wave function (79) decays in time and describes an outgoing wave for created 
particles by the Coulomb field. One may understand the quasi-stationary states by a simple 
analog of the quantum system having a potential well and finite potential barrier. In our 
case the real part of the effective potential consists of the attractive inverse square potential 
and a repulsive Coulomb barrier for qQ > 0. Therefore, the quasi-stationary state has a 
certain probability to tunnel the barrier and decays in time. Such radial wave function is 
found 

^lir) = Nr^'+^e'^'n\I{l' + 1 + iX' , 21' + 2, -2i^fer), (86) 
where M is the confluent hypergeometric function and 



A' = ^, l' = -\ + ^^^'Q'-{l + \? (87) 



Prom the asymptotic form 

g-i7r(/'+l+iA') 



-2i^/~er)-^ -'-'^ +— — ; —{-2i^'^ 



r(i' + 1 - iX') ^ ^ ^ r(i' + 1 + iX') 

(88) 

we see that, assuming Y^e{^/e) > 0, the first term corresponds to an incoming wave and the 
second term to an outgoing wave. We now require only outgoing waves and this condition 
is prescribed by 

I' + l-iX' = -rir, (89) 

which makes first term vanish due to r(— n^) = ±oo. Then we obtain the complex energy 
eigenvalues 

u; = m\l+( ^ ]] . (90) 



.n-{l + l)+i^q^Q^-{l + ly 
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We make a few comments. The consistency of quasi-stationary states requires that 
Re(a;) > 0, Im(a;) > and Re(-v^) > 0, which can be shown by a straightforward alge- 
bra. The result does not depend on the relative sign of charges, and Eq. (90) holds true also 
for the opposite charges, qQ < 0. The complex energy eigenvalues (90) are the analytical 
continuation of the real energy eigenvalues of bound states when qQ > / + |. Moreover, 
the complex energy eigenvalues are quantized and indeed describe the quasi-stationary or 
meta-stable states that decay into infinity. More importantly, the pair production rate given 
by the imaginary part of frequency is proportional to m-\/gQ. This implies that the stronger 
the electric field is, the more the pairs of charged particles are produced. 



B. Charged Black Holes 

Now the radial equation for each mode is described by Eq. (80). At sufficiently away 
from the black hole {r/M ^ 1), r* ~ r and Eq. (81) becomes approximately 

vM « ^^^JQ_211 + ^(^±11^. (91) 

Hence, Eq. (80) is approximately the same as that of the Minkowski spacetime with the 
replacement of uo by uqQ — M. Therefore, we are able to obtain the wave function valid in 
the asymptotic regions 

^;(r) = Nr^'+^e'^'^^Mil' + 1 + iXb, 21' + 2, -2i^/^r), (92) 

where 

A, = 5^?^, r=-i+.^W-(( + l^ (93) 
The outgoing wave is found by imposing the condition 

/' + 1 - = -rir, (94) 

which removes the incoming wave. Finally, after some algebra, we obtain the complex 
frequencies 



UJ — m 



l+(f) ]-'^\ MqQ 



1 - M2F2{l + (f ) }J F2{l + (f ) } 



where 



F^n- (l-r\)-r iJq'Q' - {l + . (96) 
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In the limit of M = 0, we recover the result (90) of the Minkowski spacetime. 

We make a few comments. There is a formal similarity of Eq. (80) and the Regge- Wheeler 
equation for the axial perturbations of the black hole [20, 21]. The boundary condition for 
the purely outgoing wave leads to a discrete spectrum of complex eigenfrequencies for both 
the quasi-stationary states and quasi-normal modes [24, 25, 26, 27]. However, the quasi- 
stationary states exist even in the Minkowski spacetime, whereas the quasi-normal modes 
vanish in the Minkowski spacetime. This is ascribed to the fact that the pair production 
is a universal phenomenon occurring in all spacetimes whereas the quasi-normal modes are 
a consequence of the spacetime curvature. More importantly, the decay rate, given by the 
imaginary part of the frequency, has an additional contribution in proportion to the mass of 
the hole. This implies that a charged black hole is a more efficient mechanism for Schwinger 
pair production. 

C. Pair Production 

Finally, we apply the canonical method in Sees. Ill and IV to calculate the pair- 
production rate from charged back holes. At the spatial infinity r +oo, (r* +oo), the 
efi'cctive potential vanishes, Vi{r* = oo) = 0, and at the outer horizon r — >• r_)_, (r* — > —oo), 
it has 

Vi{r+)^u;^-m^-(u;-^Y. (97) 
V r_|_ ^ 

Thus there is a potential barrier between r* — — oo and r* — +oo. We find the instanton 
action 




where Tq and are the turning points of Vi{r*). Then the mean number of produced boson 
pairs is 

M = e-'^'('^). (99) 

The mean number (99) is the leading term for the pair-production rate and was derived also 
in Ref. [28]. The pair-production rate for bosons per volume per time is given by 

2Im£^ = (2a + 1)^(2/ + l)ln(l + e-2^'('^)) 

l,UI 
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(2(7+1) 



E (-1) 



oo 



n+1 



{21 + 1) 



(100) 



2n 



;=0,n=l 



n 



and for fermions by 



2Im£{ = -(2(7 + l)5^(2; + l)ln(l-e-2^'('^)) 



27r , ^ , n y 

t=U,n=l 



(101) 



VI. CONCLUSION 

We studied the pair production by inhomogeneous electromagnetic fields, in particular, 
by charged black holes. We applied the recently developed canonical method to find the pair- 
production rate. For a static uniform electric field, the field equation can be solved either in 
the time-dependent gauge or in the space-dependent gauge. The concept of particle creation 
by an external field applies directly to the field equation in the time-dependent gauge. This 
time-dependent gauge, however, may not be practically applicable for inhomogeneous fields 
due to the mixed nature of potential. In the case of the space-dependent gauge, the mode- 
decomposed field equations have tunneling barriers, and the tunneling interpretation may be 
used for pair production. In this interpretation the tunneling probability is related with pair 
production and the no-tunneling probability with the vacuum persistence, i.e., the vacuum- 
to- vacuum transition. This canonical method can readily be applied to even inhomogeneous 
fields. Further, the instanton action including all order corrections leads to an accurate 
formula for the pair-production rate. 

Finally, we studied pair production by charged black holes. The Klein-Gordon equation 
in a charged black hole background has the form of the Regge- Wheeler equation. For a 
strong electric field, the mode-decomposed field equations have the quasi-stationary states 
that describe the purely outgoing waves for produced pairs. Their complex frequencies 
determine the decaying rate for the quasi-stationary states. We also noted a remarkable 
similarity between the quasi-stationary states of created charged particles and the quasi- 
normal modes of the axial or polar perturbations of a black hole. The same boundary 
condition for purely outgoing waves is used in both cases and, as a consequence, both wave 
functions have a discrete spectrum of complex frequencies. We then applied the canonical 
method to find the pair-production rate in terms of the instanton action by charged black 
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holes. 
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